Exact RG Flow in an Expanding Universe and Screening of the Cosmological Constant 
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We obtain the exact renormalization group (RG) flow equation for a self-interacting real scalar 
field in an expanding cosmological background. The beta functional for the potential in the local 
potential approximation is determined in terms of the mode functions defining the vacuum of the 
theory. For a scalar in an inflating universe, which is released in the Bunch-Davies vacuum, the 
flow equation is analyzed in detail using the optimized cutoff function for the quartically truncated 
potential. In that case, the constant term in the potential, which can in fact be viewed as the 
cosmological constant, decreases as the IR cutoff is removed, supporting the existence of a quantum 
mechanical screening mechanism. On the other hand, the running of the quartic coupling constant 
at IR is altered by the expansion, which seems to affect the quantum triviality of the theory. 



I. INTRODUCTION 

The exact renormalization group (ERG) method is an 
important tool that can be used to probe nonperturba- 
tive properties of quantum field theories (for review see 
e.g. [1-6]). In the form derived in [7], the ERG flow 
gives an interpolation between the UV (or fundamen- 
tal) action describing the microphysics and the full (IR) 
quantum effective action, which encodes the macroscopic 
phenomena. The method offers an effective way of calcu- 
lating the beta functions and the nontrivial fixed points 
of RG flows. The ERG is obtained by introducing an IR 
cutoff function in momentum space, which is arbitrary 
as long as it satisfies certain conditions. Although the 
physical results become independent of the cutoff func- 
tion after the regulator is removed, some specific choices 
are known to optimize the convergence and the stabil- 
ity of the flow [8, 9]. In general, the flow equation is 
an infinite-dimensional functional differential equation, 
which can be simplified in various ways and/or trunca- 
tions. For instance, in the local potential approximation 
(LPA), which is first introduced in [10, 11], one may take 
a momentum-independent limit of the ERG with con- 
stant field values. It is also possible to employ a system- 
atic derivative expansion [12] (or expansions in powers of 
fields or canonical dimensions, see e.g. [1]). The standard 
perturbation theory can be recovered from the exact flow 
by an iteration in the coupling constant [13, 14]. 

The ERG techniques are applied to pure gravity and 
gravity with matter systems in searching for nontrivial 
UV fixed points implying asymptotic safety [15]. There 
are also suggestions that the long-distance gravity is de- 
scribed by an IR fixed point [16-19] (see also [20-26] 
which discuss possible implications of the running cos- 
mological and Newton's constants). Recently, a concrete 
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cosmological model is proposed, which takes into account 
the RG running of gravitational and matter coupling con- 
stants [27] (see also [28], which takes the cutoff to be 
proportional to the Ricci scalar giving an effective RG 
improved f(R) model ). 

In this paper, we consider a self-interacting real scalar 
field propagating in a fixed cosmological background and 
obtain the ERG flow equations for the in-out and in-in 
quantum effective actions. Our method differs from the 
above-mentioned papers in that we try to determine the 
impact of the cosmic evolution on the quantum dynamics 
of the scalar field and ignore any possible backreaction ef- 
fects. In other words, our approach fills the gap between 
the flat space studies and the works that include the grav- 
itational dynamics. After introducing a cutoff function 
in the comoving momentum space that is suitable for a 
scalar in a cosmological background, the Lorentz signa- 
ture ERG equations can be obtained both for the in-out 
and in-in formalisms. Unlike the flat space case, the ERG 
equations must be interpreted with care due to the ex- 
plicit time dependence of the background. Following [29], 
we assume that self-interactions are turned on for a fi- 
nite time affecting the running of the coupling constants. 
This interpretation is viable for inflation where one is in- 
terested in the form of the effective action at the end of 
the exponential expansion. 

Although wc first derive exact flow equations using a 
comoving cutoff, it is also possible to obtain the ERG 
equations with a physical cutoff by redefining the flow 
parameter. Since the physical scale deviates from the co- 
moving scale in time, the coarse graining in the effective 
action is different for both cases. While intuitively cut- 
ting off the physical scale seems the correct procedure, 
in the loop calculations the comoving IR cutoff can be 
argued to be more natural [30] . We will see that the two 
flows give qualitatively similar results in de Sitter space. 

As usual in the ERG analysis, one needs to employ 
some form of approximation to get information from the 
functional differential equation. The simplest route is to 
apply the LPA and analyze the flow for the effective quan- 
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turn potential. In this approximation, we find that the 
beta functional for the potential is given in terms of the 
held mode functions dehning the vacuum of the theory. 
We then consider an inflating universe and truncate the 
potential with a quartic polynomial after which the flow 
equations become a system of ordinary nonlinear differ- 
ential equations for three parameters; the cosmological 
constant, the scalar mass and the quartic coupling con- 
stant. The cosmological constant term in the potential is 
usually discarded in the flat space analyses but we keep it 
to understand its running with the renormalization scale. 
The flow equations become tractable in certain limits, 
which provides a qualitative understanding of the scale 
dependence of the physical parameters. 

The plan of the paper is as follows: In the next section 
we obtain the ERG flow equations both for the in-out and 
the in-in quantum effective actions. We then apply the 
LPA and determine the beta functional of the potential 
in terms of the field mode functions. In our approxima- 
tion, the in-out and in-in cases become identical, which 
is similar to the one loop calculations of the quantum 
effective potential. In section III, we focus on inflation 
and truncate the potential up to quartic order to study 
the runnings of the relevant parameters in detail. In the 
conclusions we summarize our findings and comment on 
some open problems for future work. 



II. ERG FLOW IN A COSMOLOGICAL 
SPACETIME 

In this paper we consider a self-interacting real scalar 
field (j), propagating in a cosmological background that 
has the metric 

ds 2 = a{n) 2 {-di 1 2 +dx 2 ), (I) 

where n is the conformal time. The scalar has the stan- 
dard action: 

S[4>] =~\J d^x^—g [V^W + V(<j>)] . (2) 

In flat space one usually writes the ERG equations in 
the Euclidean signature. Since it is not possible to Wick 
rotate a general cosmological spacetime, here one needs 
to work with Lorentz signature. 

In obtaining the ERG flow, the action (2) is modified 
by employing an IR cutoff function i? M (fc 2 ), which we 
take it to depend on the comoving momentum variable k 
and a comoving renormalization scale [i. As in flat space, 



we require R IJi (k 2 ) to obey 1 

Rnik 2 ) ->• as k -> oo, 

R^k 2 ) -> p? as k -> 0, (3) 

Rpi(k 2 ) — > oo as /j, — > oo, 

R^(k 2 ) -J-O as 0, 

and introduce the scale-dependent action as 

SM = lJd V d 3 ka 2 (> 2 - [k 2 + R^k 2 )]^ 2 - a 2 v) , 

( 4 ) 

where the prime denotes an 77-derivative and <f> is the 
Fourier transformed field in the momentum space 

We will use this tilde convention for all Fourier- 
transformed variables below. Since (j> is real (p(k)* — 
4>{—k) and k dependence of <j> in (4) is suppressed. 
Let us introduce the in-out generating functional 

W-fDi^M+if*,*,*) (6) 
and the quantum effective action 

ry<fl - wyj] - J d v d 3 xj^ - as^], (7) 

where A5 M is the extra cutoff-function-dependent piece 
added to the action (4), i.e. 

AS M = -\J d V d 3 ka 2 R^k 2 )0 2 . (8) 

In general, to make the in-out path integral (6) well de- 
fined, one must specify in and out asymptotic regions 
that would determine the necessary boundary conditions. 
As an apparent manifestation of this problem, one must 
actually fix the limits of 77- integrals in (6). For the mo- 
ment we ignore this important technical point by assum- 
ing that suitable conditions are imposed to make (6) well 
defined and discuss the issue in more detail below. 

Before writing down the ERG flow equation for the 
quantum effective action T^^], let us briefly elaborate 
on the conditions (3). In flat space, each of these condi- 
tions is imposed for a reason and one should make sure 
that they have the same consequences in this cosmologi- 
cal setup. While the first condition in (3) ensures that the 
cutoff function does not change the UV behavior of the 
theory the last condition is needed to remove the cutoff 



1 As pointed out in [1], the third condition must ensure the diver- 
gence of the cutoff function i? M (A; 2 ) in such a way that a delta 
functional arises in the path integral to kill all quantum fluctua- 
tions, see the discussion below. 
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completely, i.e. r M=o [0] gives the full quantum effective 
action. 

The second condition in (3) must provide that R ll (k 2 ) 
works as an IR regulator. In flat space this is guaranteed 
by the fact that the corresponding term in the action 
actually becomes a mass term for IR modes. Although 
this is no longer true for the action (4), since the mass 
term has a different geometric structure (note the scale 
factor dependencies), it is not difficult see that R fl (k 2 ) 
functions as an IR regulator. Consider, for example, a 
massive scalar in de Sitter space. In the free theory, the 
mode function of the field is given by the Hankcl function 
4>(k) = H n (-kr))/^-wr]/A where n = ^9/4- m 2 /H 2 . 

Therefore, as ki] — > one has <j>(k) ~ 1/fc™, and for a 
light field with real n, this singular behavior may give 
strong IR effects (or divergencies for the massless field). 
When the free action is modified as in (4), the new mode 
functions can be obtained by k 2 — > k 2 + R^(k 2 ) and as 
krj — > one has 4>{k) ~ 1//U n due to the second condition 
in (3). This 2 makes the theory IR safe. 

Finally, the third condition is imposed to satisfy 
1^=00 [0] = S[4>] so that r^[</>] interpolates between the 
classical action and the full quantum action as /i runs 
from oo to 0. To verify this, one may use the functional 
integral representation of T M [<f>] , which can easily be de- 
rived from (6) and (7) to read 



= / D X exp[iS[<p + x}-i I X^f 



J d V d 3 ka 2 R^k 2 )x(k) : 



(9) 



In the limit p — > oo, the divergence of R^ (k 2 ) must ensure 
that the second line produces a delta functional 3 
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which then implies r M=00 [(/>] = S[<f>] by forcing \ — 
0. For instance, similar to the n-dimensional identity 

lim^ 0O cxp(-^ 2 X; 4 a; 2 /2)Wy2^)™ = Hi^i), i* is 
enough to have i? M (fc 2 ) — > p 2 as p — > oo. Note that 
the delta functional is defined so that in a path integral 
/ Dx^M-^Ix] — ^[O] f° r an y well-defined functional F[x] 
and the extra constant terms arising from the delta func- 
tional limit can be absorbed in the path integral measure. 
We therefore see that the IR cutoff function introduced 
in the action as in (4) has all the desired properties. 

To derive the ERG flow equation, one may differentiate 
(7) with respect to p. The derivative of W P [J] can be 



2 Note that the argument works for a comoving cutoff. If the cutoff 
is chosen to be physical then the mode equation is modified in a 
nontrivial way. 

3 To be more precise, the delta functional arises if R /J- (k 2 ) has a 
small imaginary piece. Otherwise one may argue that the integral 
becomes infinitely oscillatory which gives a nonzero contribution 
only when \ = 0. 



found from (6) as 

dWjM _ 
dp 



l - f d V d 3 ka 2 ^ <4> 2 > . (II) 
2 J dp 



On the other hand, the derivative of A5 M is 

%M = _1 [ dr jd 3 ka 2 ^<4» 2 , (12) 
dpi 2 J dfi 

where we used the fact that in the Legendre transfor- 
mation (7), the classical field <fi equals the vacuum ex- 
pectation value of the quantum field operator. One thus 
obtains 



>^--ld V d* k a 2dR » S2W ^ J] 



dT,j\<j)) _ i_ 
~ 2 



dfi 



dp 5J(ri,k)5J(ri,-k) 



, (13) 



where the source J must be viewed as a functional of 
<j>. Note that the factor of i in (13) arises since <<fi 2 > 
— <4>> 2 = -iS 2 W/SJ 2 . Since <j>, and thus J, are real, 
the Fourier-transformed variables are not independent. 
To avoid any complication, the functional derivatives in 
(13) can be defined as 



5 J{r), k) (2 



1 [ffi xe &M 



5J(r],x) : 



(14) 



which can be verified from the Fourier transformation 
and the chain rule. Therefore, (13) can be rewritten as 

dY u -i 



dp 2(2tt) 3 



/ 



drj d 3 xd 3 yd 3 k a 2 ^^-^^G^, x; rj, y), (15) 



(10) where the Green function is defined by 



*' $?$f$& **™ ^ s <"'->")^-fl. 

(16) 

Eq. (15) is an exact functional differential equation for 
the quantum effective action T M , which gives ERG flow 
in a cosmological space. Note that in principle the Green 
function is determined by T M [(/>] from (16). 

It is possible to repeat the above calculation for the in- 
in formalism. In that case the the generating functional 
is given by 

e *wyj+,j-] = J £ ) 0± e »[s M [^+]-5 ( .[^-]+/(j+^ -j- 4>-)] j 

(17) 

where Dcjr^ denotes the in-in path integral measure [31] 
and the modified actions are given by (4) both for <f> + 
and <p~ . The quantum effective action is defined by the 
Legendre transformation 

T„[4>+,4>-\ = w,[J + , J-] - J (J+4>+ - J-0-) 

-(AS^+l-AS^-]). (18) 
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There is a finite upper time limit rjf for the in-in path 
integral and the sum in (17) is over all fields with 
0+fo) = ^-(r, f ) (and 4>'+{ nf ) = #-{ Vf ), see [31]). 
Therefore SW^/SJ+^f) = -SW^/ '5J~ (r?/), which im- 
plies that the effective action T ll [(j) + , cj>~] is only defined 
for cj) + {rjf) = (f>~(t]f). Differentiating (17) with respect 
to jU and using (18) one may obtain 



dfj, 



1 = 1 
2 

5 2 W a 



dr]d 3 ka 2^E 

d\x 



(19) 



S 2 W U 



SJ+(r], k)SJ+{r], -k) 6J~{r], k)SJ~(r], -k) 

where must be viewed as functionals of ± and the 
functional derivatives in the momentum space should be 
understood as in (14). Defining 



f(r),x; r)',y) 



s 2 r^ 5 2 ry 

SJ+(r).x)SJ+( V ',y) SJ+( V ,x)SJ-( V ',y) 
SJ~ {rj,x)5J+(r}' ,y) 8J~ (t],x)SJ~ (rj f ,y) 



where 

and the matrix Green function 



/ 



dijd 3 zT(r]', x; ij, z)IG()j, z; 77", y) = I5(r/' 



where 



" 1 




'10" 


-1 


, 1 = 


1 



the flow equation (19) can be written as 



2(2tt) 3 



(20) 

■rj")6(x-y), 
(21) 

(22) 
(23) 



d V d 3 xd 3 yd 3 k a 2 e **-(*-tf ^Tr [lG(r], x; 77, y)} , 

where Tr denotes the two-dimensional matrix trace. Eq. 
(23) is the ERG flow equation for the in-in quantum ef- 
fective action since the Green function G is in principle 
fixed by T (/>-] from (21). 

As mentioned above, the asymptotic boundary condi- 
tions must be specified to define in-out and in-in path 
integrals properly. However, once the ERG equation is 
obtained, one may forget about its origin and view it as a 
property of the renormalized and finite quantum effective 
action. On the other hand, one still needs to impose cer- 
tain boundary conditions to solve in-out and in-in Green 
functions, (16) and (21), uniquely. As we will see, the 
boundary conditions are directly related to the vacuum 
chosen for the theory. 

From the first property in (3), the momentum inte- 
grals in (15) and (23) are guaranteed to converge at in- 
finity (the three-dimensional spatial integrals correspond 
to Fourier transformations and thus they are not prob- 
lematic). But r\ integrals may diverge if the asymptotic 



regions and the fall off conditions for the fields are not 
chosen carefully. To avoid any further complications one 
may follow [29] and assume that the system is released at 
time 7]i in a specific vacuum state, which can be associ- 
ated with a past infinity that is imagined to exist before 
r/i. For the in-in case this would fix a finite interval from 
this initial time rji to the final time rjf at which the in- 
in correlation functions are calculated. Similarly, for the 
in-out case one may assume that the interactions (and 
the expansion of the space) are turned off at some time 
rjf, which can then be joined to a future asymptotic re- 
gion. As a result we simply suppose a finite interval for 
the time integrals, which can be determined according to 
the physical problem at hand. For example, in applying 
the ERG for a scalar in an inflationary background, jjj 
and rjf can be taken to be the beginning and the ending 
of inflation. In that case the ERG gives how the expo- 
nential expansion affects quantum dynamics at a given 
renormalization scale. 

The quantum effective action can be systematically ex- 
panded in the number of field derivatives. For the in-out 
case, one may write 



drjd 3 xa 2 [Z((j>,fi)((j>' 2 - dc/> 2 ) - a 2 V(<p,n)] 

(24) 



where the dotted terms represent fields with at least four 
derivatives. Note that the general covariance fixes the 
form of the two derivative terms, i.e. the wave-function 
renormalization Z(4>,n) is single. On the other hand, 
Z(cj),fj,) and the effective potential V(<j),n) are expected 
to depend on the conformal time explicitly since the back- 
ground is time dependent. The kinetic operator can be 
calculated from (24) as 



6 2 r L 



54>{x)5<t>{y) 



2a' d 



-a 2 Z{^n)[— + ——-d 



where 



drj 2 ' a drj 



x=(r],x), y = (v',y), 



(25) 
5(x - y), 

(26) 



and the dots in (25) denote the terms containing at least 
two derivatives of the field <f>. Till now no approxima- 
tion is used and to proceed one must calculate the Green 
function satisfying (16). Obviously, the inverse of (25) is 
hard, if not impossible, to calculate for arbitrary <f>. In- 
deed, even the operator that is first order in the derivative 
expansion, i.e. the second line of (25), is difficult to in- 
vert 4 for a generic <j). However, it is possible to evaluate 



4 To invert (25) in the derivative expansion one may use the for- 
mula (A + B)- 1 = A^ 1 - A~ 1 BA~ 1 + ... where A represents the 
second line of (25) and B corresponds to the higher-derivative 
terms indicated by the three dots. 
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(15) at a constant field. In that case, to calculate (16) 
it is enough to consider the second line of (25), which is 
very similar to the free kinetic operator other than the 
possible explicit conformal time dependencies of Z(<j>, fi) 
and V{4>, fi). 

To find the Green function (16), which is evaluated 
at a constant field 0, one may construct a free quantum 
operator as 



fa) 



i 



ay/Z{(j>,n) 



d 3 k 



(27) 



where H.c. means Hermitian conjugate and the mode 
functions obey 



k2 + R,(k 2 ) _a^ + 



d 2 V 



Z(<j),n) a 2Z(<f>,fj,) d(j> 2 



Xk = 0. (28) 



The free quantum operator <j)(x) should not be confused 
with the constant field <$>, which is used in the evaluation 
of (15). With the help of 0, the Green function can be 
constructed as 



G(x,y) = -i<T(4>(x)<l>(y)) >, 



where T denotes time ordering. Although the field opera- 
tor <f) is introduced to calculate the Green function, there 
is a close connection with <f> and our original quantum 
scalar field 4>: In perturbation theory, the two fields ac- 
tually coincide with each other at the lowest order. This 
allows us to fix the arbitrariness in the Green function, 
which shows up as the arbitrariness of the solutions of 
(28). Namely, in calculating (29), one should use the so- 
lution of (28) consistent with the vacuum chosen for the 
theory. Keeping this in mind, (15) reduces to 



/ 



drj a 



A dV(<j>,n) 



d\x 



1 

8^ 



, ,i , dR u 
drjd 3 k ^ 



dfi Z{4>,n) 



(30) 



which can be viewed as the beta functional of the quan- 
tum effective potential. Let us remind the reader that 
in the flat space ERG calculations, a similar reasoning 
is used to find the exact propagator in momentum space 
where one simply takes the inverse of a polynomial ob- 
tained from the kinetic term of the effective action, which 
indeed refers to the vacuum. Since (30) is precisely (15) 
evaluated for a constant field, it is exact and free of any 
approximations. 

For the in-in case, the zeroth-order classical UV action 
is given by T[<p + ,(f>~] = S[<fi + ] — S[<j)~]. Although quan- 
tum corrections are expected to generate ± mixings in 
the effective action, we assume that 

r^+,^-] = r^+}-r^-}, (3i) 

where the functional has the form (24). The main rea- 
son for this approximation is that the exact in-in effective 



action vanishes for 



i.e. iy</>,0] = 0. Thus, the 



mixing terms can be expected to be suppressed (indeed 



in perturbation theory the mixing does not arise at one 
loop). As a crucial technical point, if one considers a 
general form for r M [0 + , cf>~], then uniquely defining the 
matrix Green function (21) becomes a problem since one 
may no longer refer to the free theory. On the contrary, 
with the assumed form (31) the operator T defined above 
(20) becomes similar to the free in-in kinetic operator and 
one may calculate the matrix Green function G by using 
the standard in-in formalism. In that case it is easy to see 
that the beta functional for the in-in effective potential, 
which is equal for + and — branches, becomes identical 
to (30). Note that although one uses the time-ordered 
and the anti time-ordered products in defining the diag- 
onal entries of the in-in matrix Green function, this does 
not make any difference in the flow equation (23) since 
the Green functions are evaluated at equal times. 

To proceed, we conveniently choose the optimized cut- 
off function [8, 9] 

R^(k 2 ) = Z{</>, v)(p 2 - k 2 )6{n 2 - k 2 ). (32) 



It is easy to see that (32) satisfies all conditions in 
(3). The presence of the wave function renormalization 
(29) Z(cf>, /i) in (32) is required to preserve the scaling behav- 
ior of the scalar field (see e.g. [1]). Using (32) in (28) 
one sees that the mode equation becomes independent of 
k for the regularized IR modes 



k < n ■ x'L + 



d 2 V 



a 2Z{(j>,n) d(j> 2 



X» = 0. (33) 



Thus the momentum integral in (30) can be carried out 
to yield 



where rj^ is the anomalous dimension 

d\nZ{<j>,n) 



V4> = -- 



(35) 



Recall that in solving (33) one should refer to the vacuum 
chosen for the theory. 



The Flow with a Physical Cutoff 

The cutoff parameter fi refers to a comoving energy 
scale. If one would like to regulate physical IR modes 
then one should change n with time. In that case a new 
physical energy scale a can be introduced so that 



= a(rj) a. 



(36) 



The cutoff function R fl (k 2 ) depends on a through (36) 
and it also changes with time. In this new setup, a can 
be used as the new flow parameter. 
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It is possible to repeat the derivation of the beta func- 
tional, this time taking the a derivative of the quantum 
effective action in (13). This gives 



/ 



drj a 



,dV(<f>,a) _ 1 



da 



8tt 3 



J drjd 3 k 



dR» \xk{y)\ 2 

da Z(4>,a) 



(37) 



which trivially replaces (30). In this new scheme the 
optimized cutoff function becomes 

R a (k 2 ) = Z(<f>, a)(a 2 a 2 - k 2 )6(a 2 a 2 - k 2 ) (38) 

and using (38) in (28) gives the new mode equation. Un- 
like (33), the mode equation becomes piecewise defined 
in time for given constants a and k: 



a > 
a < 



k 
a 
k 
a 



X.c 



x'L + 



2 2 a " fl2 

aa -~ + 2ZU, 



d 2 V 



, 2 a 

k 2 + 



a 



a) d(j) 2 
d 2 V 



Xo 



a 2Z{4>,a) 



X* = 0.(39) 



This equation must be solved in two different regions 
and the pieces should be matched at a = k/a. The mode 
function must be continuos and should have a continuos 
first time derivative. Even though the equation in the 
range a > k/a is free from k, the matching generates k 
dependence. As a result it is not possible to carry out 
the momentum integral in (37) without solving for Xa 
explicitly. Note that the anomalous dimension rj^ is the 
same for the flow parameters \x and a. 

One may think that the flow equations with comoving 
and physical cutoffs should agree with each other since 
they look like the reparametrization of the same curve. 
This is not correct because they correspond to coarse 
graining of different scales in time and the system does 
not have time translation invariance. The exact quantum 
effective actions must be the same when the cutoffs are 
removed but the flows might differ when an approxima- 
tion is employed. In flat space, the comoving cutoff is 
analogous to a time-dependent scheme which eliminates 
different physical scales at different times. In the next 
section, we analyze both of the flow equations in detail. 
As we will see, while the comoving cutoff scheme faces 
the trans-Plankian problem, the physical cutoff scheme 
has a peculiar IR behavior. 



III. ERG FLOW EQUATIONS DURING 
INFLATION 

In this section, we analyze (34) and (37) for a scalar in 
an inflationary universe and take 



a = — 



Hrj' 



(40) 



where H is the Hubble constant. The inflation is assumed 
to last for N e- folds and typically TV > 60. We apply 



the LPA, 5 namely we take V(<j>, fi) to be independent of 
r\ (recall it may have an explicit time dependence) and 
ignore the possible contributions of anomalous dimension 
by setting 



Z{(f>,fj) 



V4> = °. 



(41) 



which leads to a closed differential equation system in- 
volving only the parameters of the potential. The LPA is 
known to be a good approximation for many applications 
in flat space [1]. 

We start with the comoving cutoff scheme. In the 
LPA, (33) becomes the free-field mode equation with a 
modified mass term and the comoving momentum value. 
Therefore we take 



X»=Xk D (k^v), 



(42) 



where xf D 1S the mode function in the Bunch-Davies 
vacuum. This gives 



X l > = \l J £ L e i H n+ VHW{-m). (43) 



where 



n = 



19 1 d 2 V 

4 ~ m 2 !^ 



(44) 



and H^p is the first Hankel function. Without loss of any 
generality one may take r\i = —1/H (so that a(%) = 1) 
and write 



-N 



(45) 



where, as noted above, N denotes the number of c-folds. 
The complex conjugate of \n can be found as 



Xi, 



I^ e - if (n + I )jff (2) ( _ w)) (46) 



which is valid both for real and complex n. 

Using (43) in (34) and carrying out fy-integrals yield 



dV 



H[i 2 



[A n (f,/H)-A n ^e- N /H)], (47) 



dn 8ir(e 3N - 1) 

where the function A n is given by 

A n (x) =x 2 [J n (x) 2 + Y n (x) 2 

-J„_i(x)J„ +1 (.t) - F„_!(a;)F„ + i(a;)] . (48) 

Here J and Y are the Bessel functions of first and second 
kinds, respectively. A n (x) is an increasing function of its 



5 In flat space, the LPA amounts to take the potential indepen- 
dent of four-momentum and to neglect the anomalous dimension. 
Since we are not Fourier transforming the time direction, in our 
case the LPA is equivalent to ^-independence of the potential. 
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FIG. 1. The graph of A„{x) for n = 0, 1/2, 1,3/2, from top 
to bottom, respectively. For n < 1, the function diverges as 
x -» 0. 



argument for real n, with a somehow weak dependence 
on the index n (see Fig. 1). For later use let us note the 
following asymptotic behavior: 



, . . Ax 1 - An 2 
A n {x) = — + — 

TT ZTTX 



as x — > oo, and 

An cot(n7r) 



A n (x) 



9 - AOn 2 + 16n 4 
32-kx 3 



0(x) 



O 



1 



(49) 
(50) 



_2» r 4 "[rw 2 -r(n-i)r(n + i)] 2 

+2 I X 



0(x 3 



as x — y for Re(n) > 0. 

Denning a new dimensionless parameter as 

—N A* 

k = e — , 

and assuming > 1, the flow equation becomes 



dn ~ 8n K [An[ 



(51) 



(52) 



which is the beta functional of the potential, i.e. one may 
Taylor expand both sides of (52) in <f> to read the beta 
functions. Note that the right hand side of (52) depends 
on the potential through the index n given in (44). 

Before analyzing the solutions of (52) for a quartic po- 
tential, we would like to comment on a few important 
points. Recall that \i has been introduced as a comoving 
scale, r/i is chosen such that a(i]i) — 1 and a (77/) = e . 
Thus, u corresponds to the physical scale in the begin- 
ning of inflation and k is the physical scale at the end of 
inflation. Since physics after inflation is our main con- 
cern, k is the scale of interest and UV/IR limits must 
refer to the variable k. 

In a realistic scenario N > 60 or so, and (52) clearly 
exhibits the trans-Plankian problem [32-34]. Namely, 
when k is near Planck scale, there appears another scale 
in the problem e N k which is trans-Plankian. On the 



other hand, for k near IR, e N k can be UV and thus (52) 
also shows UV/IR mixing, which is a peculiar property 
of de Sitter physics [29, 35]. 

As noted above A n (x) is an increasing function of its 
argument. From the basic construction of the ERG flow, 
one must imagine starting the flow at some UV scale and 
ending up at an IR scale (in principle at k = 0). We 
therefore see from (52) that V decreases 6 along the flow 
for any given value of (j). 

It is crucial to emphasize that the quantum action is 
defined once and for all in the interval {i]i,r]f) and the ac- 
tion itself does not evolve in time. Therefore, in the above 
construction the ERG flow should not be interpreted as 
an evolution in time. However, due to the external time 
dependence coming from the background, the parameters 
in the effective action are expected to vary with time. 

It is interesting to consider an almost massless field so 
that m(n) 2 /H 2 = d 2 V/dcf) 2 \^ =0 / H 2 ~ 0, which actually 
requires an extreme fine tuning of the parameters. In 
that case, one may set n = 3/2 and integrate (52) for 
6 = to obtain 



A(«) 



H 1 



[e N (K 4 ~4)-2(K 2 -K 2 )}+A , (53) 



where A = U|0 = o and Ao is an integration constant cho- 
sen such that A(k ) = A . Thus the ERG flow gives 



A(0) ~ A - f^ffV 



8tt 2 



0) 



(54) 



which shows that the the bare cosmological constant Ao 
is screened by a huge negative amount. 

The above result shows how ERG might be used to 
address some crucial details in the cosmological constant 
problem. Imagine that one calculates quantum correc- 
tions to the cosmological constant coming from a field 
theory in an expanding universe. In such a calculation 
the easiest way to regulate infinities is to introduce a 
cutoff, which is somehow ambiguous since the scales red- 
shift with the expansion. In our case, for example, in- 
flation generates a huge hierarchy of scales in time and 
the physics may depend sensitively on how the cutoff is 
chosen. The ERG equation (52) solves this ambiguity, 
since it actually gives a cumulative effect overtime and 
as a result e N factor in (54) arises. 

Let us now truncate the potential as 



V = A 



m 2 (j) 2 



A0 4 . (55) 

To get the beta functions for A, m 2 and A, one must ex- 
pand (52) in powers of (f> and equate the coefficients. In 
flat space studies, the cosmological constant term A is not 
considered in the ERG flows. In the hierarchy of parame- 
ters, A is the lowest order member, i.e. its flow is affected 
by other parameters but it does not affect the flow of any 



This can already be observed in (34) for tja < 5. 
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parameter. Therefore, ignoring it in flat space causes no 
problem other than the fact that the bare cosmological 
constant must be fine tuned to neglect any backreaction 
effects. Here, aiming to apply the ERG techniques to 
the cosmological constant problem, we keep A to see its 
scale dependent running. Since we ignore backreaction 
effects, an implicit fine tuning is also assumed so that A 
becomes small compared to the background cosmological 
constant. 



UV Regime 

It is difficult to analyze the beta functions arising from 
the flow equation (52) exactly, mainly because of flow 
varying indices of the Bessel functions. However, it is 
possible to simplify the equations in certain limits, which 
helps to understand the qualitative flow behavior. In 
the UV regime for n » 1 (recall that we have already 
assumed e N >• 1) using the asymptotic formula given in 
(49) one finds 



K > 1 



H 4 e 



e 4 _ _ 2 2, 1 

~2^ K -^2 m K + 1^2 



rn 



_ 9 3H o 3 9 \ 
Kd K m ~ --— t Xk + —^m A, 



(56) 



Let us note that in obtaining (56) no condition on to 2 
is imposed, i.e. it can be any real number including the 
negatives. From (51), the H dependence in the above 
equations becomes redundant. The flow equations for 
to 2 and A arc similar to the flat space counterparts and 
the only imprint of the expansion of the universe is the 
term that appears in the beta function of A. 
The last equation in (56) can be solved as 



•M K ) — 



A 



l-9^1n(«/«o) 



which is identical to the flat space behavior. The usual 
"triviality" argument follows from the fact that (57) im- 
plies A(/t) — > as k — > 0. However, in our case (57) 
is valid for k ~> 1 and as we will see shortly the flow 
becomes quite different at IR. 

Since A runs logarithmically for k 3> 1, one may take 
A ~ A in the second equation in (56) to get 



2 3Ao 



H 2 k 2 



I- 

K 



3A /4n- 2 



2 i 3A 2 
m + ^2 H 



(58) 

where to 2 (ko) = m 2 ,. This is again consistent with the 
flat space behavior, i.e. the renormalized mass quadrati- 
cally depends on the UV scale. Note that as one lowers 
k, m 2 gets larger. 

From (58), one finds that m 2 ~ 0(k 2 ) and using this 
in the first equation in (56), the first term in the right 



hand side can be seen to give the largest contribution. 
Integration then yields 

A(«)~|J[e J V-Kg)]+Ao, (59) 

which is consistent with (53) for k >• 1. Therefore, in the 
UV regime the greatest contribution to the cosmological 
constant is independent of the bare parameters in the 
potential, it is in some sense purely geometrical and arises 
due to the expansion of the universe. 



IR regime for a light scalar 

Consider now the next interval in the flow where one 
has 1 3> k 3> e~ N . When N is large, this interval in- 
cludes even deep IR excitations. For generic values of 
to 2 , the expansion of the function A n in (52) involves the 
derivatives of the Bessel functions with respect to their 
arguments and the expressions become very complicated. 
However, the formulas simplify considerably for a light 
scalar, which has \m 2 {n)\/ H 2 <C 1 since the index n can 
be expanded around n = 3/2. 

In that case, using (50) one finds 



nd K \ ~ [oi + a 2 ln(«;) + a 3 ln(/-c) 2 ] A 2 



(60) 



where the constants involve factors like the Euler's num- 
ber coming from the expansion of the Bessel functions 
and they are given by oi ~ 40, 02 — 21 and a 3 = 24. 
The flow of A changes drastically compared to (56), and 
(60) shows that 



A — > const, as k — > 



(61) 



invalidating the triviality argument. Indeed, the varia- 
tion of A is small since k< 1 in (50) and one may take 



(57) A(/t) <~ const. 



A 



1 + ln(«o) 



l> K >e- iv , (62) 



where A is matched with (57) at k = 1. 

In the same interval, the ERG flow equations for to 2 
and A become 

Kd K m 2 ~U 2 K 2 \b x + b 2 ln( K )]A 

+K 2 [b 3 + b A ln(K) + 6 5 ln(K) 2 ]TO 2 A (63) 



and 



nd K A ~— \e N K 4 + k 2 } + H 2 n 2 [ Cl + c 2 ln( K )]m 2 

+k 2 [c 3 + c 4 1ii(k) + c 5 ln(«;) 2 ]TO 4 , (64) 

where hi and Ci are numbers of order unity (note the ex- 
pansion in powers of m 2 /H 2 ). In the beta function of 
A, we see that the first term containing e N initially gives 
the largest contribution, but as k is lowered k 2 terms be- 
come dominant. Since the right-hand sides of (63) and 
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(64) are small in Hubble units, there appears no signif- 
icant contribution to m 2 and A in this IR region with 
1 > k > e~ N . 

Finally one may consider the extreme IR regime with 
K -C e~ N . Using the asymptotic formula (50) one sees 
that the beta functions for m 2 and A become identical to 
(60) and (63), where all the terms in the right hand sides 
are multiplied by —1. Similarly for A, one gets (64) with 
no e N term and all others in the right hand side change 
sign. These can be seen by observing that while in (52) 
A n (e N K) is much larger than A n {n) for 1 » k » e _JV , 
the opposite is true for k -C e"^. As a result, the flow 
becomes trivial and the runnings are irrelevant in the 
deep IR regime with k -C e~ N . 



Physical Cutoff 



60 




0.5 



1.0 



1.5 



2.0 



2.5 



FIG. 2. The graph of B n {x) for n = 0, 1/2, 1,3/2, from left 
to right, respectively. 



To analyze the beta functional in the physical cutoff, 
one must first solve (39). For convenience, we scale a as 



Ha 



(65) 



so that it becomes dimensionless. The solution of (39) 
can be found as 



x /^e^ + l)^(-M:.<-f (66) 

c^-kiq) 1 / 2 ^ + c 2 {-kn) l / 2 + r : r) > -f 



where n is given by (44) and 



1 d 2 V 
2H 2 deb 2 ' 



(67) 



The normalization of the solution for r\ < —a/k is fixed 
by referring to the Bunch-Davies vacuum. Demanding 
that the mode function and its first derivative to be con- 
tinuous at r\ = —a/k fixes the integration constants as 



Cl 



C-2 



7T 

16fc 
16fc 



e i§(«+5)a/ 



-<«(«) + #«(«)] (68) 
'-H'W(a)+HW(a)] (69) 



where prime denotes derivative with respect to the argu- 
ment. 

Using (66) in (37) and carrying out the momentum 
integral, one sees that the n integral in the right-hand side 
of (37) becomes J dr/a 4 , which is exactly the same with 
the left hand side. As a result, the time integrals on both 
sides cancel out each other and (37) becomes independent 
of rji, r]f, and hence the number of e-folds N. The ERG 
equation with the physical cutoff does not suffer from the 
large hierarchy of scales during inflation. The modes are 
eliminated gradually in time and this solves the trans- 
Plankian problem. 

After some algebra we find that (37) becomes 



dV 
da 



H 4 
48^' 



1 d 2 V 



B n (a), (70) 



where 

B n (a) = 2a 2 H^ 1 (a)H^} 1 (a) (71) 

-a(3 + 3n) (h & (a)H& («) + H^(a)H^\(a)) 

+ (9 + 6n - 2a 2 )H^\a)H^ 2 \a). 

Eq. (70) is the beta functional of the potential in the 
physical cutoff scheme. Here, B n (x) is a monotonically 
decreasing function of its argument for real n (see Fig. 
2), which has the following asymptotic expansions: 



B n {x) 



26 -8n 2 27+ 16n 2 (n 2 - 7) 

TTX 47T2 3 

3(441 - 1996?i 2 + 944n 4 - 64n 6 ) 



(72) 



647TX 5 



as x — > oo and 



B n ( X )=[ 2[8n2 - 9]c0t ^ + 0( X ) 
nil 

3[3 + 2n]4 n I>] : 



(73) 



0(x 2 



as x — > for Re(n) > 0. 

Even though the trans-Plankian problem is solved, 
(70) may suffer from peculiar IR behavior. If a 2 + m 2 (a) 
vanishes 7 at some point along the flow, then (70) may 
become singular depending on the behavior of the func- 
tion a 2 + m 2 (a) around its zero. For instance, as a — > 
one finds that 



dV 
da 



C- 



t 4-2y<9/4~m 2 /H 2 

a 2 + m 2 1 H 2 



(74) 



where C is some constant. If m 2 vanishes faster than or 
equal to a 2 as a — > 0, i.e. m 2 ~ a v with p > 2, then 



7 Recall that one should expand both sides of (70) around <fr = 
and m 2 = d 2 V/d(j> 2 \ 4>=0 /2. 
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(74) becomes logarithmically divergent at a = 0. This 
is related to the well-known IR divergence of a massless 
scalar in de Sitter space, which shows up as a singularity 
in the RG flow. 

It is again interesting to consider an almost massless 
field before discussing the generic case below. After ne- 
glecting mass and setting <f> = 0, which gives n = 3/2, 
(70) can be integrated to yield 



A(a) 



"24tt 2 
+A , 



[(a 4 - a 4 ) + 4(a 2 - ag) + 18In(a/a + e)] 

(75) 



where A = V|^=o, Ao is an integration constant and e 
m 2 /H 2 ~ 0. Thus the ERG flow gives 



A(0) ~ A 



ff 4 
24tt 2 



[a 4 +4 a 2-181n(e)] . 



(76) 



As in (54), we see a large screening of the cosmological 
constant, of the order of H a o.q, where a ~> 1. However, 
contrary to (54) the screening becomes independent of 
the number of e-folds, i.e. it can be thought to arise 
as a result of formulating the theory in the eternal de 
Sitter space. For a massless field, e = and the flow 
diverges logarithmically, which is the IR problem men- 
tioned above. 



UV regime 

Once more the beta functional (70) is complicated to 
solve exactly and as before we first analyze it in the UV 
regime. Truncating the potential as in (55) and using 
the expansion formula (72), for a 3> 1 (and for any real 
value of m 2 ) (70) gives 



a > 1 : 
ad a A ~ 



El 

6^ 2< 



2 2 i 1 4 

Tin a + — n m 



12tt 2 



16tt 2 



H 2 3 
ad a m 2 ~ -— -Xa 2 + ——m 2 \, 
2tt z 4tt z 



(77) 



After a constant rescaling a — » \/3a, (77) becomes iden- 
tical to (56) except the first term in the beta function of 
A. Thus, in the UV regime the comoving and the phys- 
ical cutoff schemes agree on the runnings of m 2 and A. 
They also agree on m 2 and A dependent contributions to 
A-flow but they disagree 8 on the "geometrical" term. For 
«> 1 the solutions are given by (57) and (58), where k 
is replaced by ^/3a, therefore as one lowers the cutoff, A 
decreases and m 2 increases for arbitrary initial data. 



8 We checked very carefully that the first terms in the beta func- 
tions of A in (56) and (77) are different. 



IR regime for a light scalar 

In the IR regime for a -C 1, we again focus on a light 
scalar obeying \m 2 \/H 2 < 1, so that the index n can be 
expanded around n = 3/2. Using the asymptotic formula 
(73) we find from (70) that 



2 \2 



act, A : 



a 2 X 



6tt 2 (a 2 + m 2 /iJ 2 ) 3 



[162+...] (78) 



where the dots denote the terms like a 2 , a 4 , a 2 ln(a) 
etc. that vanish as a — > 0. Contrary to the previ- 
ous flow given in (60), the IR behavior of A depends 
sensitively on the running of m 2 . This is because of 
(a 2 + d 2 V/(2H 2 d(j> 2 ))- 1 term in (70) which does not con- 
tain any additive constants. For comparison note that 
in (52), m 2 only appears through the index n and conse- 
quently it is always combined with the Hubble's constant 
H in the form 9/4 — m 2 /H 2 . Therefore, (78) must be ex- 
amined with the running of m 2 , which can be determined 
from (70) as 



ad a m 2 ~ 



H 2 



a 2 X 



-162 



36tt 2 {a 2 + m 2 /H 2 ) 2 

2/2 

2a 2 | (di 



(79) 



m 
H 2 



m 
H 2 



c?2 ln(a) -I- c?3 ln(a)^ 



where d\ ~ 34, — —56 and d% — 36. 

Examining (78) and (79), one may see that there is a 
solution for which A and m 2 approach nonzero constants: 



A/ 



m 



Mr 



as a 



0. 



(80) 



Such a solution exists since the right-hand sides of the dif- 
ferential equations (78) and (79) vanish as a — > if m 2 
approaches a nonzero constant. The measurable param- 
eters m 2 R and A_r are determined by the initial conditions 
that must be supplied for the flow. This IR behavior is 
consistent with the one obtained in the comoving cutoff 
scheme. 

However, there is another interesting solution of (78) 
and (79), which has completely different IR characteris- 
tics. Assuming that 



2„2 



(78) gives 



A : 



4tt^ 
27 



1 + — ^ 

+ H 2 



a 



Using this in (79) fixes the value of m 2 , as 
ml = --H 2 . 



(81) 



(82) 



(83) 



Therefore (78) and (79) admit an IR scaling solution for 
a <C 1, which is approximately given by 



A^-a 



1 



2„2 



as a 



0. 



(84) 
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Note that in deriving (78) and (79) we have assumed 
\m 2 \/H 2 <C 1 and the scaling solution satisfies this re- 
quirement. 

On the other hand, the beta function of the cosmolog- 
ical constant A can be determined from (70) as 



ad a A ~ 



H 4 



216tt 2 a 2 + m 2 



[162+ 



(85) 



H 2 



[ei + e 2 ln(a)] 



m 



e 3 + e 4 ln(a) + e$ hi(a) 2 ] 



where ei ~ -97, e 2 ~ 108, e 3 ~ 84, e 4 ~ -102 and 
e 5 = 36. For the first solution (80), the ERG flow (85) 
can be integrated to yield 



A 



H 4 



216tt 2 m 2 R /H 2 



An. 



(86) 



As a — > 0, the change in A is negligible and the flow 
becomes irrelevant at IR. However, the integration of (85) 
for the scaling solution (84) gives 



A~ 



H 4 



1627T 1 



■ln(a) + A , 



(87) 



which shows that A receives a large (divergent) negative 
IR contribution. Since A diverges as a — > 0, one may 
prefer to keep it small but finite, which implies a phys- 
ically observable nonzero value for m 2 . If m 2 obs denotes 
this measurable value, the scaling solution (84) gives a 
minimum for a as 



a* = 2 



m' 



bs 



H 2 ' 



The change in A in the IR regime is given by 



AA ~ 



H 4 



162tt ; 



■ ln(a*), 



(88) 



(89) 



which should be added to the amount of UV screening. 
Note that when m 2 < the potential minima becomes 
degenerate at <j> min = m 2 /(2A) and there emerges 

a global Z 2 symmetry, which is broken by the nonzero 
scalar vacuum expectation value. 



IV. DISCUSSION 

In this section, we would like to elaborate on the inter- 
pretation and possible implications of our findings. As 
noted previously, to obtain the full quantum effective ac- 
tion the cutoff must be sent to zero, and thus the comov- 
ing and the physical cutoff schemes should agree with 
each other in the exact quantum picture. However, the 
action is usually truncated as we did above and conse- 
quently one is not dealing with the full quantum action 
but instead an effective description relevant to a given 
physical situation. In that case, the ERG flow parameter 
must be set to the scale of interest in the problem. 



It is important to emphasize that we take a fixed grav- 
itational background and neglect any possible backreac- 
tion effects. The background geometry sets a scale, which 
is the Hubble constant H, and in fixing the value of the 
ERG flow parameter, this scale must be taken into ac- 
count. 

In the light of the above comments, let us, for exam- 
ple, elucidate (53) and (54), which give the ERG flow of a 
massless scalar in an inflationary spacetime. First, recall 
that although k is defined to give the physical scale at the 
end of inflation, it is actually a comoving scale measured 
in units of H. Therefore (53) should be thought to give 
the effective value of A at the given comoving scale k, 
provided that A is the effective vacuum energy density 
at the comoving (UV) scale kq. This standard interpre- 
tation becomes subtle in the comoving cutoff scheme due 
to the time dependence of the background, which shows 
up as V-dependence in (53). Thus, the RG flow fixes the 
effective value of A at a given comoving scale at a given 
time and (53) is expressed to yield the scale dependence 
of A at the end of inflation. Specifically, (54) gives the 
total screening of A from UV to IR regimes at the end of 
inflation (since k is defined to give the physical scale at 
the end of inflation) and if, for instance, one would like 
to determine the impact of the vacuum energy density 
at cosmological scales, it is A and not An that must be 
taken into account. Similarly, if one is interested in the 
UV physics at scale kq at the end of inflation, the value 
of A must be used as the corresponding vacuum energy 
density. 

The RG flow (75) in the physical cutoff scheme is much 
more easy to interpret. As before, (75) gives the effective 
value of A at the physical scale a but this time A(a) turns 
out to be time independent. Here, the time dependence 
is absorbed in the flow parameter, which is chosen to be 
the physical scale redshifting with expansion. One may 
imagine that the effective vacuum energy density of the 
quantum scalar decreases as it flows from sub-Hubble UV 
scale to the super-Hubble IR scale due to the constant 
spacetime curvature acting as a potential barrier. 

The above results describe how the effective vacuum 
energy density of a test scalar changes from UV to IR 
in a fixed inflationary background and it is encouraging 
to see that the RG flow implies a fall off, i.e. the classi- 
cal gravitational field screens the vacuum energy density 
at large scales during inflation. However, at this point 
it is not possible to claim that this provides a realistic 
screening to give a graceful exit mechanism for inflation 
as in, e.g., [36]. To verify this, backreaction effects must 
be taken into account, i.e. the gravitational dynamics 
should be determined in the presence of the matter/scalar 
(and possibly metric) fluctuations (alternatively, one may 
also include the gravitational degrees of freedom in the 
ERG flow analysis to obtain an effective action including 
the gravitational degrees of freedom as in [27, 28]). 

Keeping these in mind, one may try to understand pos- 
sible implications of our findings for the cosmological con- 
stant problem. In our analysis, the RG scale is considered 
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to be a free parameter which can be set according to the 
physical problem at hand. In determining the cosmologi- 
cal impact of A, it is natural to set the physical RG scale 
at the instantaneous Hubble expansion rate as in [36, 37] 
(see also [38], which uses the scale factor of the universe 
to fix the RG scale). This motivates the physical cutoff 
scheme and let us consider the flow (75). In writing (75), 
we had scaled the flow parameter a in (65) to make it 
a dimcnsionlcss parameter. To avoid any complication, 
one may simply rescale a — > a/H to obtain 



A(a) 



1 



[(a 4 -a 4 )+4H 2 (a 2 -a 2 ) 



24tt 2 

+181n(a/a + e)]+A . 



(90) 



Note that this rescaling corresponds to a simple 
reparametrization of the same RG trajectory that does 
not change physics. Let us now assume that this equa- 
tion is still valid when H is changing (at least slowly) 
and the cosmological impact of the vacuum energy is ef- 
fective at the Hubble scale so that a ~ H. Assuming 
further that q « ao, the effective cosmological constant 
at the Hubble radius becomes 



A(H) ~ A 



24tt 2 



H 2 *l 
6ir 2 ' 



(91) 



where we ignore the logarithmic term to keep the dis- 
cussion simple. The UV scale can naturally be identified 
with the Planck scale, i.e. Ao ~ Mp and ao — Mp. The 
first subtracted term in (91) is the flat space contribu- 
tion and the second one is related to the expansion of 
the universe. As usual, the characteristic quartic quan- 
tum contribution to the cosmological constant appears 
with the opposite sign in the RG flow analysis, but the 
sign of the last geometric term is nontrivial. In principle, 
the bare cosmological constant Ao can (mostly) be can- 
celed out by the af, term by fine tuning. On the other 
hand, (91) shows that for large H the screening of A is 
large and as H decreases A increases. Since A(H) is sup- 
posed to act like a cosmological constant in Einstein's 
equations, (91) implies the existence of an "equilibrium" 
value for H (if the last term in (91) would appear with a 
positive sign, then this would indicate an instability, i.e 
growing H would give an increasing A(H), which would 
then force H to grow more via Einstein's equations) . Al- 
though the RG flow indicates a screening of the vacuum 
energy density for a massless scalar reducing the value 
of the bare cosmological constant, extra ingredients are 
necessary to solve the graceful exit problem. 

One possible modification is to consider a massive field, 
which is supposed to yield m 4 and m 2 H 2 terms in the 
ERG equations. Indeed, these terms can already be seen 
to arise in the flow (64), however they are not expected 
to give significant contributions since we have assumed 
m <C H. On the other hand, our results show that the 
ERG analysis can provide a theoretical basis for the phe- 
nomenological models studied in [39-42], which mainly 
consider cosmological consequences of a vacuum energy 
density fixed by a polynomial of the Hubble parameter. 



For small H, the screening of the bare cosmological 
constant is small by (91). This is somehow expected 
since our analysis uncovers the impact of the spacetime 
curvature on the ERG flow. Therefore, at least for a 
massless field, our results are not useful to understand 
the accelerated expansion of the universe today (to avoid 
any confusion let us note that the extreme IR regime dis- 
cussed in section III does not refer to our universe today, 
but it arises in the inflationary regime due to the expo- 
nential hierarchy of scales related to the trans-Plankian 
problem). However, it would be interesting to consider 
the derivative expansion in an expanding universe along 
the lines of [12] to determine the form of the effective po- 
tential near a fixed point, which would give a nontrivial 
dependence on the Hubble parameter. 



V. CONCLUSIONS 

In this paper, we obtain the ERG flow equations for the 
in-out and in-in quantum effective actions correspond- 
ing to a self-interacting real scalar field in an expand- 
ing cosmological spacetime, which are given by (15) and 
(23), respectively. As in flat space, the ERG equations 
can be derived using the path integral representation, 
however, formulating the theory in a cosmological back- 
ground gives rise to some key differences. Although it is a 
minor technical detail, the derivation must be carried out 
in Lorentz signature since it is not known how to Wick 
rotate a general cosmological spacetime. Moreover, one 
can only Fourier transform the three spatial dimensions. 
The ERG equations involve the exact propagator that 
should be obtained from the quantum effective action by 
(16) and (21). In flat space, at least in some approx- 
imations, this is a trivial procedure where one simply 
takes the inverse of a polynomial in momentum space. 
In the expanding universe, the Green function must be 
constructed by referring to a vacuum state. 

May be the most striking difference compared to the 
flat space ERG calculations is the need for an extra in- 
put for the interpretation of the quantum effective ac- 
tion. One usually thinks of an action to be defined for 
the whole spacetime. Moreover the fields must be im- 
posed to obey certain boundary/fall off conditions for 
the variational principle to work. In the absence of nice 
asymptotic regions as in a generic cosmological space- 
time, defining a well behaved action and a variational 
principle can be problematic. To avoid technical com- 
plications, we follow [29] and assume that the cosmo- 
logical background is only active for a finite time interval 
{r)i,r)f), which is glued to suitable past and future asymp- 
totic regions. Supposing further that the interactions are 
turned off outside this interval, the action functionals, 
including the quantum effective action, can be thought 
to de defined in the region (jjj,??/). Such a description 
is suitable for inflation, where r\i and r\f mark the be- 
ginning and the ending of the exponential expansion. In 
that case, the quantum effective action encodes the im- 
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prints of the classical gravitational field on the quantum 
dynamics of the scalar. 

Yet another difference with the flat space ERG analysis 
is the existence of comoving and physical cutoff schemes. 
In ERG theory, the term scheme usually refers to a choice 
of a cutoff function. In an expanding universe one can 
consider two main coarse graining procedures with re- 
spect to the comoving or the physical scales. Although, 
full quantum effective actions must agree as the cutoffs 
are removed, using the comoving or the physical cutoff 
schemes can give different flow behavior when an approx- 
imation is employed. In the case of inflation, we observe 
that the flow obtained with a comoving cutoff function 
exhibits the trans-Plankian problem. On the other hand, 
the IR behavior of the flow with the physical cutoff func- 
tion turns out to be richer. 

It is possible to evaluate the ERG equation, which is 
an infinite dimensional functional differential equation, 
for a constant field to obtain a beta functional for the 
quantum effective potential, which is a system of infi- 
nite dimensional ordinary differential equations. For in- 
flation, (52) and (70) give the beta functionals of a self- 
interacting scalar field, in the comoving and the physical 
cutoff schemes, respectively. In a time dependent back- 
ground, one would expect the parameters in the quantum 
effective action also to depend on time. In applying the 
LPA, we neglect any plausible time dependencies. Nev- 
ertheless, since the whole setup depends on a given time 
interval (r]i,r]f), one would expect to see its trace in the 
equations. For inflation, the interval is fixed by the num- 
ber of e- folds N , which shows up in the beta functional 
(52). One may think that the time dependent parameters 
of the exact quantum potential are replaced in a way by 
the time averaged counterparts in the LPA. The situa- 
tion in the physical cutoff scheme is more interesting. As 
noted above (70), in that case the time integrals in the 
left- and right-hand sides of the quantum effective action 
have the same form leading to (70). This shows that (70) 
is valid for any interval and one can indeed set r\i = — oo 
and r]f = covering the whole Poincare patch of the de 
Sitter space (note that (52) is ill defined for N — > oo). 
Therefore, the ERG flow with the physical cutoff can be 
used to probe the quantum dynamics of an interacting 
scalar in the eternal de Sitter space. 

As one of our main results in the present work, we 
quantitatively determine the ERG flow of the cosmologi- 
cal constant, which is viewed as the constant term in the 
scalar potential. Starting the flow at some UV scale from 
a given value of the cosmological constant, we find that 
the it's value decreases as the cutoff is lowered. Not sur- 
prisingly, the largest decrease occurs in the UV regime, 
which is of the order of CH A where C is a very large num- 
ber depending on the scheme and the UV cutoff, which 
can presumably be taken as the Planck scale. For the 
comoving cutoff, the running of the cosmological con- 
stant becomes irrelevant at IR, which is also the case 
for the trivial IR solution of the physical cutoff scheme 



(80). However, there is a special scaling solution to the 
flow equations in the physical cutoff scheme in which the 
cosmological constant can be screened by a large neg- 
ative amount even at IR. These findings are consistent 
with the expectations that the cosmological constant is 
relaxed by strong IR quantum effects [43-46] (or similarly 
by the backreaction of long wavelength cosmological per- 
turbations [47-49]). Although the backreaction effects 
are neglected in this work, the emergence of screening 
in a nonperturbative calculation supports the relaxation 
mechanism. 

As another interesting finding, we observe that the 
beta function for the quartic coupling constant changes 
drastically at IR compared to the flat space behavior. 
As a result, the coupling constant does not vanish as 
the IR cutoff is removed, which invalidates the trivial- 
ity argument. This is a curious IR effect, which emerges 
due to the presence of spacetime curvature. However, the 
story is more interesting since in the scaling solution (84) 
the coupling constant vanishes as the cutoff is removed, 
although at a different rate compared to the flat space 
flow. This indicates the existence of a nontrivial phase 
diagram for the quartically interacting scalar field in de 
Sitter space. 

It is possible to extend the present work in various dif- 
ferent directions. It would be interesting to understand 
the backreaction effects along the lines of [27] or [28] to 
figure out how the screening mechanism affects the gravi- 
tational dynamics and the other way around. In addition 
to the LPA, which is occasionally insufficient in capturing 
the correct quantum behavior for some cases [1], there are 
several other approximation methods developed to ana- 
lyze the ERG flow equations in flat space, some of which 
can be naturally adopted in the cosmological framework. 
For example, considering O(N) scalar models in the large 
N approximation or applying the systematic derivative 
expansion as in [12] are appealing possibilities. Finally, 
in this paper we approximate the in-in quantum effec- 
tive action in the form (31), however it is clear that this 
approximation is not sufficient to capture all the crucial 
effects arising from the peculiar properties of the in-in for- 
malism. Besides, in determining the backreaction effects, 
one should obtain in-in quantum effective field equations, 
which can be derived from the in-in action by applying 
a (f> + functional derivative and then setting (f> + — (f>~ . It 
is an important open problem to generalize the ERG for- 
malism for a nonperturbative analysis of the in-in quan- 
tum effective field equations. 
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